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Abstract 

We present a pedagogical study of the commonly employed Speed-Plot (SP) and Time-delay (TD) 
methods for extracting the resonance parameters from the data of two particle coupled-channels 
reactions. Within several exactly solvable models, it is found that these two methods find poles 
on different Riemann sheets and are not always valid. We then develop an analytic continuation 
method for extracting nucleon resonances within a dynamical coupled-channel formulation of ttN 
and reactions. The main focus of this paper is on resolving the complications due to the 
coupling with the unstable ttA, pN, aN channels which decay into ttttN states. By using the 
results from the considered exactly solvable models, explicit numerical procedures are presented 
and verified. As a first application of the developed analytic continuation method, we present the 
nucleon resonances in some partial waves extracted within a recently developed coupled-channels 
model of vrA^ reactions. The results from this realistic vrA^ model, which includes vrA^, rjN, vrA, 
pN, and aN channels, also show that the simple pole parametrization of the resonant propagator 
using the poles extracted from SP and TD methods works poorly. 
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I. INTRODUCTION 



The excited baryon and meson states couple strongly with the continuum states. Thus 
they are identified with the resonance states in hadron reactions. The spectra and decay 
widths of the hadron resonances reveal the role of the confinement and chiral symmetry 
of QCD in the non-perturbative region. Therefore, the extraction of the basic resonance 
parameters from reaction data is one of the important tasks in hadron physics. Ideally, it 
should involve the following steps: 

1. Perform complete measurements of all independent observables of the reactions con- 
sidered. For example, forpseudo- scalar meson photoproduction reactions one needs 
to measure 8 observables [ij : differential cross sections, three single polarizations E, 
T and P, and four double polarizations G, H, E, and F. 

2. Extract the partial- wave amplitudes (PWA) from the data. Here we need to solve a 
non-trivial practical problem since all observables are bi-linear combinations of PWA; 
i.e. a ~ fus'fls- 

3. Extract the resonance parameters from the extracted PWA. Here the often employed 
methods are based on the Breit-Wigner form[2], Speed-Plot method of Hoehlerjs], 0], 
and Time-Delay method of WignerjH, 0]. A more sophisticated and rigorous method is 
to use the dispersion relations, K-matrix, and dynamical model to analytically continue 
the PWA to the complex energy plane on which the resonance poles and residues are 
determined. Extensive works based on these three models are reviewed in Ref. 0]. 

In reality, we do not have complete measurements for practically all meson-nucleon reactions. 
Even if the measurements are complete, the step 2 requires some model assumptions to solve 
the inverse bi-linear problem in extracting PWA. This model dependence must be taken into 
account in interpreting the extracted resonance parameters. 

In this work, we focus on the step 3 in conjunction with the recent efforts in extracting 
the nucleon resonances from very extensive and high quality data of meson production 
reactions, as reviewed in Ref. The nucleon resonances A^*, A listed by the Particle Data 
Group [s^ are mainly from the analysis of 'kN scattering and pion photoproduction reactions. 
The Speed-Plot and Time-Delay methods are most often used in these analyses since they 
only require the PWA determined in the step 2. The purpose of this work is to examine 
the extent to which these two methods are valid and to develop an analytic continuation 
method within a recently developed dynamical model[^ of meson production reactions in 
the nucleon resonance region. 

It is useful to first briefiy recall how the resonances are defined in textbooks. By an- 
alytic continuation, the scattering T-matrix can be defined on the complex energy plane. 



Its analytical structure is well studied [10|, lUl, 112, ll3|, for the non-relativistic two-body 



scattering. For the single-channel case, the scattering T-matrix is a single-valued function of 
momentum p on the complex momentum p-plane, but is a double- valued function of energy 
E on the complex energy E-plane because of the quadratic relation p = |2mi?|^/^e"^^/^. 
Therefore the complex E-plane is composed of two Riemann sheets. The physical (p) sheet 
is defined by specifying the range of phase < 0^; < 27r, and the un-physical (u) sheet by 
2vr < < 47r. As illustrated in Fig. [T|, the shaded area with Imp > of the upper part of 
Fig. 1(a) corresponds to the physical i?-sheet shown in Fig. 1(b). Similarly, the unphysical 
ii^-sheet shown in Fig. 1(c) corresponds to the Imp < area in the lower part of Fig. 1(a). 
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FIG. 1: The complex momentum p-plane ((a)) and its corresponding complex energy E'-plane 
which has a physical sheet ((b)) and unphysical sheet ((c)). Their correspondence is indicated by 
the same number. Solid squares (circles) represent the bound state (resonance) poles. 



On the physical sheet, the only possible singularities are on the real E-axis : the bound 
state poles (solid square in Fig. 1(b)) below the threshold energy E^h and the unitarity cut 
from Efh to infinity. On the unphysical sheet, a pole (solid circle in Fig. 1(c)) on the lower 
half plane and Re Epoie > Eth corresponds to a resonance. From unitarity and analyticity of 
S-matrix, each resonance pole has an accompanied pole, called conjugate pole, which exists 
on the upper half of the un-physical sheet. A resonance pole is due to the mechanism : an 
unstable system is formed and decay subsequently during the collision. The mathematical 
details of this interpretation can be found in textbooks, such as in chapter 8 of Goldberger 
and Watson[15i]. 

For multi-channel case, the analytic structure of the scattering T-matrix becomes more 



complex [111, [12, ll3|, [IJ]. We postpone the discussion on this until section III where a 
two-channel Breit-Wigner form of scattering amplitude will be used to give a pedagogical 
explanation. 

The essential point of the Time-Delay and Speed-Plot methods is that the resonance 
poles discussed above can be determined from the partial-wave amplitudes defined on the 
physical real energies. The concept of time delay was originally introduced by Eisenbudjsl] 
and discussed by Wignerj^, Dalitz and Moorehouse[3l, and Nussenzveig[l3]- It is defined 
as the difference between the time in which a wave packet passes through an interaction 
region and the time spent by a free wave packet passing though the same distance. It was 
generalized to the time delay matrix in multi-channel problems and discussed further in 
terms of a lifetime matrix by Smith [l8|. Kelkar et al.|l9l. 20| applied this multi-channel 
formulation to develop a practical Time-Delay method to extract hadron resonances. 

The Speed-Plot method was developed by Hoehler0, 15] to extract the nucleon resonances 
from the vrA^ partial-wave amplitudes. It is also based on the concept of time delay, but he 
identified resonances by "speed" which is the absolute value of the energy derivative of the 
scattering amplitude. The speed is always positive by its definition while the time delay can 
be negative and becomes "time advance" [2l|. 
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It is generally assumed that both the Speed-Plot and Time-Delay methods give approxi- 
mate positions of the resonance poles which are rigorously defined by the analytic continu- 
ation of scattering amplitudes. However their validity is not clear at all. Moreover it is not 
clear which Riemann sheet the pole positions extracted from using these methods belong 
to. In this paper we analyze several exactly solvable models to clarify these questions. In 
particular, we find that these two methods give poles on different Riemann sheets. 

We then develop an analytic continuation method for extracting the resonance parameters 
within a recently developed Hamiltonian formulationj^ of meson-nucleon reactions. We first 
establish the method by using several exactly solvable models. The main task here is to 
handle the singularities associated with the unstable particle channels such as ttA, pN , and 
aN . We then apply the method to extract the nucleon resonances in some partial waves 
of ttN scattering within the dynamical coupled-channels model developed in Ref.j22|. The 
analytical structure of the resonance propagator is also analyzed within this model. 

In section II, we describe the formula for applying the Speed-Plot (SP) and Time-Delay 
(TD) methods. These two methods are then analyzed and tested in section III by using 
the commonly used two-channel Breit-Wigner form of S-matrix. In section IV, we develop 
an analytic continuation method by using several exactly solvable resonance models and 
further test the SP and TD methods. Section V is devoted to resolve the complications 
due to couplings with unstable particle channels. In section VI, we present the results from 
applying the developed analytic continuation method to extract the nucleon resonances in 
some partial waves within the vriV model of Ref.H- A summary is given in section VII. 



II. FORMULA FOR TIME-DELAY AND SPEED-PLOT METHODS 

The objective of the Time Delay (TD) and Speed-Plot (SP) methods is to determine the 
mass Mr and width Tr of a resonance from the S-matrix of reactions. In this section, we 
will not explain how these two methods were introduced, as briefiy discussed in section I. 
Rather we only give their formula in practical applications. 

For the single- channel elastic scattering case, the time delay of the outgoing wave packet 



with respect to the non- interacting wave packet in a given partial wave is defined [6|, Il6|, [17 



as 



where the S-matrix is related to the phase shift 5 by 

S{E) = e^'^ . (2) 
Eqs. (P) and lead to a simple expression of time delay 

At = 2-. (3) 

The time-delay (TD) method is to find the resonance mass M^j by finding the maximum 
of the time delay 

dAt , 

-dE \'-''-= 
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As an ideal example, we evaluate the time delay for the S-matrix defined by the well known 
Breit-Wigner resonance formula 



Eq. ([T]) then leads to 



S{E) 



E-Mr- iTr/2 
E-Mr + iTr/2'- 



(5) 



At 



R 



{E-MRf + VllA' 



(6) 



which obviously takes a maximum aX E = Mr and hence Mr is defined as the resonance 
mass. It also gives the following simple physical interpretation of the width Tr in terms of 
the time delay of the wave packet passing through the interaction region 



At I 



E=Mr 
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(7) 



Here and in the rest of this paper, the normalization is chosen such that the S-matrix is 
related to the T-matrix by 



S{E) = 1 + 2iT{E) 



(8) 



We then note that for the S-matrix Eq. ([5]), the width can also be expressed in terms of 
T-matrix as 



2 



\T\ 

\dE\ 



(9) 



E=Mii 



In the analysis of vrA^ scattering, Hoehlerjl, 0] introduced a speed plot (SP) method. The 
speed is defined as 



sp{E) 



dT 



dE 



Using Eqs. (El)-©, Eq. leads to 



sp(E) 



lAtl 



(10) 



(11) 



Thus the speed is also related to the time delay of the wave packet. The SP method defines 
the resonance mass Mr by the maximum of the speed 



^Sp(E) \E=Mn-- 



0. 



(12) 



From Eqs. (jl]), ffTTl) . and f|T2l) . it is obvious that the speed plot and time delay will give the 
same resonance mass for the single-channel elastic scattering case. 

For the multi-channel reactions with N open channels, the S-matrix becomes a. N x N 
matrix S{E). Smith 18f| introduced a life-time matrix Q defined by 



Q{E) = -tS{E) 



^dS{E) 



dE 



(13) 



5 



The above equation can be considered as an extension of Eq. ([T]) of the single- channel 
case. The trace of the life-time matrix can be expressed in terms of eigen phases 5i of the 
S-matrixjii, H 111] 



TtQ{E) = 2- 



dE 



The resonance mass is then obtained by finding 

dTrQ{E) 



dE 



\e=Mr-- 



0. 



(14) 



(15) 



However, the eigen phases Si can be obtained only when we know all of the S-matrix elements 
associated with all open channels. In practice, only the elastic scattering amplitude and a 
few of the inelastic amplitudes can be extracted from the data. Therefore Eqs. (fT^ - (fT5|) of 
Smith's TD method can not be used rigorously in practice. 

Here we focus on the TD method used by Kelkar et al. 1^, 2^, 2l[. They defined the time 
delay for the channel i only by the diagonal component 5*^ of the S-matrix and its derivative 



td{E) 



Re -i 



1 dSii 
Sii dE 



(16) 



Obviously, this method is identical to Eq. ([T]) of the single channel case except that the 
S-matrix element here is Sa = rje^"^^ with rj denoting the inelasticity. Eq. f|T6|) can be 
considered as an approximation o f Eg. (fT3l) b v neglecting the inelastic channels in summing 
the intermediate states. In Refs. 
maximum of the time delay 



|20| . |21| , the resonance mass Mr is determined by the 



dE 



td{E) 



\E=Mr 



and the width Tr by 



td{Mn±TR/2) = td(Mfi)/2. 



(17) 
(18) 

(19) 



The SP method for multi-channel case is simply to define the speed by the diagonal Ta 
of the T-matrix 



sp(E) 



dE 



(20) 



and define [sj] the width by assuming that the T-matrix element can be parametrized as 

Tr/2 



Tu{E) = n{E) + z{E)- 



(21) 



E-Mr + zTr/2' 

Here Th is an non-resonant amplitude and the resonant amplitude is defined by the reso- 
nance mass Mr, the width Tr, and a complex residue z{W). By assuming that the energy 
dependence of Th(W), z{W), and Tr can be neglected at energies near Mr, Eqs. ( !20!) and 
(!2T|) obviously satisfies Eq. (fT2l) and lead to the following condition 

sp{Mr±Tr/2) = sp(Mr)/2. (22) 

Eqs. f l20l) . f|T2l) and fl22l) are used in applying the SP method to extract the resonance mass 
Mr and width Tr from the partial-wave amplitudes. Eq. fl2T|) is not needed in practice, but 
is an essential assumption of SP method. 
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III. ANALYSIS OF SPEED-PLOT AND TIME-DELAY METHODS 



To examine the SP and TD methods, we consider a commonly used two-channel Breit- 
Wigner (BW) amplitude which can be derived ll|, |26|, |27| from the analytical property 
of the S-matrix. To make the contact with what we will discuss in the rest of this paper, 
we will indicate here how this amplitude can be derived from a Hamiltonian formulation of 
meson-baryon reactions, such as that developed in Ref.[^. 

It is sufficient to consider the simplest two-channel case with a non-relativistic two-particle 
Hamiltonian defined by 

H = Ho + V. (23) 

In the center of mass frame Hq can be written as 

2 

Ho = Yl K) Ki + + (24) 

where rriik is the mass of k-th particle in channel i, and /Xj = miimi2/{mii + mi2) is the 
reduced mass. In each partial-wave, the S-matrix is a 2 x 2 matrix and can be written 

^ = TT^' ^^^^ 

where p is the density of state, and the K-matrix, which is also a 2 x 2 matrix, is defined by 
the following Lippmann-Schwinger equation 

K{E) = V + V^^K{E). (26) 

Here P means taking the principal-value of the integration over the propagator. 

We now consider the on-shell matrix element of the S-matrix Eq. (J25l) . If the on-shell 
momentum is denoted as pi for channel i, we then have 

± m pK\j) = 6ij ±iiTpiKij, (27) 

where pi = Pifii- The S-matrix element of the 1 — 1 elastic scattering is then of the following 
explicit form 

^ _ (1 - i7ipiKii){l + i7ip2K22) - T^"^ P1P2K12K21 

(1 + mpiKii){l + mp2K22) - -r"^ P1P2K11K22 

If we assume that at energies near the resonance energy the K-matrix can be approximated 

as 

K . ^ _9iE^ /29) 
where M is a mass parameter which is a real number, Eq. ( l28l) can then be written as 

= (30) 

E - M + zpi7i + ip2l2 
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Here we have defined 7^ = TigffJ.i > 0. If we furtlier assume that 7i is independent of 
scattering energy, Eq. (l30l) is the commonly used two-channel Breit-Wigner formula j26l. [27I . 



30[ |. In the rest of this section, we will follow these earlier works and treat 71 and 72 as 
energy independent parameters of the model. 

Since the scattering T-matrix is related to the S-matrix by 

Su{E) = l + 2tTn{E), (31) 

Eq.(|30|) leads to 

T{E) = T^,{E) 

-^'P' (32) 



E - M + i'jipi + i^2P2 
From nowon we use the notation T[E) for the 1 — >^ 1 amplitude Tii{E). 



A. Analytic Properties of the S-matrix 

Within the two-channels Breit-Wigner model specified above, we will analyze in this 
subsection the analytic properties of the S-matrix on the complex energy E-plane. This will 
also allow us to explain clearly some terminologies which are commonly seen but often not 
explicitly explained in the literatures on resonance extractions. 

The on-shell momenta pi for channel i is defined by 

where 

= (mil +mi2). (34) 
We can define the threshold variable A between two channels by 

E,= P^ + ^, (35) 



where 



2/ii 



77221 + 1^22 - mn - mi2 (36) 



is the threshold energy of the second channel. 

The momenta at poles of the S-matrix Eq. fl5U|l can be determined by solving 

E-M + ipi-ii + ip2i2 = 0. (37) 

By using Eqs. the above equation can be written as 

pt + apl + hpl + cpi + ci = 0, (38) 
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where 



a = iAfiiji, (39) 

6 = 4/ii(^272 -^-mi), (40) 

c = -SifilMji, (41) 

d = A{ixlM'-ix,fi2liA'). (42) 



Eq. fl38l) means that the Breit-Wigner amphtude Eq. fl30|) has four poles. Each pole is 
specified by two on-shell momenta Pa = {pia,P2a) with a = 1,2, 3, 4. The analytic properties 
of the S-matrix Eq. fl5Ul) depends on how these poles are located on the complex energy 
iJ-plane. As we explained in section I, the energy plane for each channel has two Riemann 
sheets because of the quadratic relation Eq. (1351) between the momentum pi and energy Ef, 

namely Pi = ^2/Xj|i?j|e*'^*/^ for i = 1,2. For each channel, the physical (p) sheet is defined by 
specifying the range of phase < (pi < 2n, and the un-physical (m) sheet by 2tt < (pi < in. 
The correspondence between the momentum pj-plane and the energy ii^j-plane is similar to 
that illustrated in Fig. [TJ For the considered two-channel case, we thus have four energy 
sheets specified by the signs of Impi and Imp2: PP, up, uu, and pu, as shown in Fig. 2. Thus 
each of four poles Pa = {pia,P2a) can be on one of these E-sheets. 

To be more specific, we now consider the case which is most relevant to our study of 
nucleon resonances. That is the KeEi > and Re i?2 > case that the poles are all above 
the thresholds of both channels. From Eq. fl37j) . we immediately notice that if {pia,P2a) 
with E = Ea is one of the solutions, {—pla: ~P2a) with E = E* is also a solution. Therefore 
the four poles determined by Eq. fl38l) can be grouped into two pairs. In the following 
discussions, they are denoted as {Ea,El) and {Eh,E^). Without losing generality, one can 
assume that one of the poles is in the range of (Repi > 0, Rep2 > 0) and the other in the 
range of (Repi < 0, Rep2 < 0). If the first pole {pia,P2a) is in the region where (Repi^ > 0, 
RGP2a > 0) and {Impia < 0, Imp2a < 0) , it is a pole, denoted as Eji, on the uu-sheet 
of Fig. [21 This pole is usually called the resonance pole and is closer than other poles 
on up or pu-sheets to the physical pp-sheet, as will be explained later. In the Hamiltonian 
formulation considered in this work and the well developed collision theory, a resonance pole 
can be mathematically derived p^] from the mechanism that an unstable system is formed 
and decay subsequently during the collision. The resonance pole Er has an accompanied 
pole E"^ at {—play ~P2a) which is also on the -u-u-sheet as shown in Fig. [2l E^ is called the 
'conjugate pole' of Ej^. 

The second pole at {pib,P2b) with Impi;, < and Imp2f) > ( Impn, > and Imp26 < 0) 
may be on the wp-sheet (pu-sheet), depending on the parameters 71 and 72. This pole is 
called the shadow pole|28|. A shadow pole on up-sheet and its conjugate pole are Es and 
E*s in Fig. H 

We now note that in this simple BW model, the zeros of the S-matrix Eq. fl30l) . where 
Sii{E) = 0, is defined by its numerator 

E-M - ipai + ip272 = 0. (43) 

The above equation can be cast into the form of Eq. (137]) by simply replacing pi by —pi. 
Thus solutions of Eq. (1431) . called the zeros of S-matrix, can be readily obtained from the 
solutions {pia,P2a) and {pib,P2b) of Eq. ([37]). They are {-pia,P2a) with Ea and (p^^ -pD 
with i?* for a = a,b. The zero at {—pib,P2b) is on the pp-sheet, denoted as Ezs and E^s in 
Fig. [2j Similarly, the zero at (— pia,P2a) is on the pw-sheet, shown as Ezr together with its 
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Im E 
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Ezs X Re E 





Re E 
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(a) (b) 



Im E 
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Im E 
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Es O Re E 





Re E 








EzR X 


up-sheet 




pu-sheet 





(c) (d) 

FIG. 2: Poles and zeros of the simplified two-channel Breit-Wigner form (fii = ^2 and A = 0) 
of S-matrix on the complex E-plane which has pp , uu , up and pu sheets. The open circles on 
the uu-sheet ((b)) are the resonance pole Er and its conjugate pole i?^. The open circles on the 
up-sheet ((c)) are the shadow pole Es and its conjugate pole E"^. The crossed on the pp-sheet ((a)) 
are the zero Ezs a-nd its conjugate E'^g which are at the same energies of the shadow poles Es 
and Es- The crossed on the pu-sheet ((d)) are the zero Ezr and its conjugate -E^/j which are at 
the same energies of the resonance poles Er and E"^. 



conjugate ii^^^ in Fig. [2J Note that Fig. [2] is for the case that the parameters 71 and 72 are 
chosen such that the shadow poles Es and its conjugate are on the up-sheet. For other 
possible 7i and 72, the pole positions could be different from what are shown in Fig. [21 but 
their close relations, as discussed above, are the same. 

From the above analysis, it is clear that the poles and zeros of the S-matrix are closely 
related. Their locations on the 4 Riemann sheets can be conveniently displayed on one 
complex plane by introducing a variable t 2^ 



Pi = (44) 

- 2a/|^, (45, 

Hence each point in the t— plane corresponds to a set of (^1,^2)- In Fig. [3l the reso- 
nance position Er and shadow Es poles and their conjugate poles and zeros of S-matrix 
Ezs, Ezr, E*^g, E^pi^ are shown on t-plane. The physical S-matrix at real energies, which 
determine the observables, are on the bold lines. The zero energy and the threshold of the 
second channel correspond to t = z and t = 0, respectively. One can see that the resonance 
pole Er is closer than the shadow pole Es to the bold lines (S-matrix) and hence can have 
the largest effect on the observables. Consequently, most of the rapid energy dependence of 
observables are attributed to the resonance poles, not the shadow poles or the other poles 
shown in FigjSl On the other hand, the zero E'^^ of the S-matrix is also close to the bold 
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FIG. 3: The poles and zeros of S-matrix shown in Fig. 2 are displayed on the complex t-plane 
defined by Eqs. and (pD . 



lines. As seen in the derivations given above, this zero E^g is closely related to shadow pole 
Es- Thus the shadow poles can also be related to the observables. Of course, which pole is 
most important in determining the rapid energy dependence of observables also depends on 
the residues of the T-matrix at the pole positions. 

In the next subsection, we will use a further simplified BW form to explain more clearly 
the close relations between the poles and zeros of the S-matrix. In particular we will see 
explicitly that the shadow pole Es in Fig. [2], which is not on the same Riemann sheet as 
En, can be located by the zeros of the S-matrix. More importantly, we will also see how the 
SP and TD methods work both analytically and numerically. 

B. Positions of poles 

For simplicity, we assume that the threshold energies of the two channels are the same 
and hence ^ = fii = fi2 and A = 0. Therefore we have simple relations between energy 
and momenta : E = = pl/2fi and p = pi = ±P2- As discussed in the previous 

subsection, the four Riemann sheets are classified by the sign of the imaginary part of the 
momentum; namely, physical (unphysical) sheet is assigned by Imp > (Imp < 0). With 
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the simplification p = pi = ±p2, we obviously have pi = p2 = p on the pp and -uw-sheets, 
and pi = —p2 = p on up and pu-sheets. 

Let us start with the case of pi = p2 = p for the pp-sheet or -u-u-sheet. The S-matrix 
element Eq. (1301) of the first channel can be written as 

It can be cast into the following more transparent form 

SniE) = + (47) 

ip-PR)iP + PR) 

with 



Pn = J2fiM - {fi^+y - ifij+, (48) 



PS = V2/iM - (/x7_)2 - 2;i7_, (49) 

where 

7± = 7i ± 72. (50) 

Remembering that we consider 71 > and 72 > 0. To make use of Fig 12] in the following 
discussion, we consider the case that 71 > 72 and hence 7± > and both p^ and ps defined 
in Eqs. (I48|) - (l49|) are associated with the unphysical w-sheet. For the case of 7_ < 0, is 
associated with the physical p-sheet and the following presentation can be easily modified 
to account for this case. 

Clearly, Eq. (147|) means that the S-matrix has a pole a.t pi = p2 = p = Pr on the uu-sheet 
with a resonance energy 

2 

En = ^ = M-ft^l- 27+v//^(2M-^7i)- (51) 

Its conjugate pole E'^ is at pi = p2 = p = —p*ii- The positions of Er and E*^ are shown 
in the upper right side of FigO Eq. flTTl) also indicates that the zero of S-matrix is at 
Pi= P2 = —Ps which is on the pp-sheet because of Im(— pg) > 0. The energy of this zero of 
S-matrix is 

2 

Ezs = Y^ = M-^^^l- Z7_ {iI{2M^^^). (52) 

Its conjugate E^g is at pi = P2 = p*s- The positions of Ezs and E^g are on the pp-sheet, as 
shown in the upper left side of FigI21 

We next consider the pi = — p2 = p case that the poles and zeros of the S-matrix are on 
the up or pw-sheets. The S-matrix Eq. (1301) for this case then takes the following form 

> ^ _ M + ^(7l - 72)p ^ ' 
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By comparing Eq. (H6l) and Eq. (l53|l and using the variables pr and ps defined by Eqs. ( HHj) 
and (H9|) . Eq. (153|1 can be written as 



The above equation indicates that for the considered 7_ > 0, the S-matrix has a shadow 
pole at pi = —p2 = p = ps on the wp-sheet. Thus its position Es = p|/(2/i) is identical to 
Ezs of the zero of the S-matrix on the pp-sheet; namely 

Es = Ezs 

2 

= = M - /i7! - 27-\/M2M-/x7^). (55) 

This means that the shadow pole Es on the wp-sheet can be found from searching for the 
zero Ezs of S-matrix on the pp-sheet. 

Eq. fl5^ also gives a zero of S-matrix at pi = —p2 = p = —pn on pw-plane because 
Im(— p/j) > 0. Its energy Ezr is also identical to Er defined above 

EzR = Er 
2 

= 1^ = M - /i7^ - 27+V/^(2M-/i7i). (56) 

The positions of Es and Ezr and their conjugates i?^ and E^r are also in the lower parts 
of FigH 

From the above analysis for the 7_ > case, we see that the energies of the resonance 
poles may be obtained by studying the poles of the S-matrix on the uu-sheet and those of 
the shadow poles may be obtained from the zeros of the S-matrix on the pp-sheet. The 
analysis for the 7„ < case is similar. Here we only mention that when 7_ > is changed 
to 7_ < 0, the shadow poles Es and Eg on the up-sheet move to the pu-sheet and zeros Ezs 
and E'^g will be on the ww-sheet. 

Now let us examine how the SP and TD methods can find the poles defined by the above 
exact expressions of the two-channel BW amplitude. We first recall that in applying the 
SP and TD methods, the energy E and momentum pj in the S-matrix are restricted on the 
positive real-axis. For the considered simplified case, we thus have pi = P2 = P and the 
T-matrix Eq. (!32|) then becomes 

T(E) = . (57) 

^ ' E-M + i-i+p ^ ' 

Our task is to examine whether the resonance mass Mr and width Y r found by applying 
the SP and TD methods on the expression Eq. flSTl) are close to the real and imaginary parts 
of the poles defined in the previous subsection. 

According to Eqs. f|T2|) and fl20|) . the SP method finds the resonance mass Mr by finding 
the maximum of the speed through the use of the condition 

dE dE „ 

L J E=M[i 
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With the T-matrix Eg. (1571) . Eq. (l58i) leads to the following equation 



3M| + 



(3M + 2/i7^)M^ - (7M^ - 6fxM-fl)MR + = 0. 



This equation can be written in the dimension-less form as 

3E^ + (3 + 2a)E^ - (7 - 6a)E + 1 



0, 



(59) 



(60) 



with E = Mfi/M,a = fijI/M^. With some inspections, one can see that Eq. fl^U]) has 
real and positive Mn solutions only in the a < 0.417 region. Two of the three solutions are 
the maximum and minimum points of the speed, and the third one is less than 0. The SP 
method defines the maximum point of the speed as the resonance mass. We find that this 
solution can be expanded as 



R SP 



M 



1 



AM 



2 



M2 



(61) 



Clearly, [MrIsp equals to the real part of Eq. ( 15T1) if we neglect the second and higher order 
terms in the expansion in powers of h'jI/M. Therefore the SP method is accurate only 
under the condition that fi'yI/M ^ 1. Moreover, it is clear from the above equation that 
speed has no stationary point for /i7^/M > 0.417 and therefore the SP method will fail to 
find the pole even if there is a pole within the model. 

We next turn to discussing the width obtained by the SP method. It is evaluated by 
using Eq. IQ. We find that it can also be expanded as 



[r 



RlSP 



\T\ 




dT 




dE 


[Mr]sp 



/x7i(2M-/x72) 



M 



(62) 



Here again, if we neglect higher order terms of fi'-fl/M, the SP method can give the imaginary 
part of Er of the exact expression Eq. (I5T1) . As we have seen in Figj2l the two-channel BW 
S-matrix has two pairs of poles on the uu-sheet and the wp-sheet. However the SP method 
can only find the pole Er on the -uw-sheet. 

From the above analysis, it is clear that the accuracy of SP method is controlled by 
/i7^/M. We examine this by using an example with /i = J^'^^ , frtN = 938.5 MeV and 
= 139.6 MeV, and M = m^v + + 600 MeV. In FigJlfthe solid curves are the pole 
positions on the ww-sheet. They are obtained from evaluating the exact analytical formula 
Eq. flSTl) for 7-/7+ = 0.5 and varying yU7^/M from 0.01 to 0.6. With the same parameters 
we then apply the SP method to search for pole {Mr, —iTR/2) from the amplitude Eq. (|57|) 
numerically by using Eqs. ([9]), ffTOl) and f|T2|) . The obtained pole positions were the filled 
squares shown in FigJH As expected, the SP method works very well for small /i7^/M. 
However the pole position from SP starts to deviate from the exact results (solid curves) as 
/i7^/M increases. There is no filled squares in FigJUin the region fij^/M > 0.417 because 
the SP can not find a pole in this region. This is not because of the numerical accuracy of 
our calculation, but is the intrinsic limitation of the SP method as, discussed above. 

We now turn to investigating the TD method. We apply Eq. f[T7|) to search for the 
resonance mass Mr from the amplitude Eq. (!57|) . Since it is not clear how to interpret 
Kelkar's prescription Eq. (fT9|) . we assume that the S-matrix element Su is of the following 
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FIG. 4: The dependence of the pole positions on the un-sheet (soHd curve) and up-sheet 

(dashed curve) of the simpUfied two-channels Breit-Wigner form (^1 =^2)^ = 0) calculated at 
7_/7+ = 0.5. The real (imaginary) parts of poles are shown in (a) ( (b)) of the figure. The solid 
squares ( solid circles) are the results obtained from using the SP (TD) method. Triangles in (b) 
are the widths of the poles obtained using Eq. p9]) . 



form 



Sii{E) 



E-M - iTs/2 
E-M + iTR/2 



Then Eq. ( fT9l) leads to an improved expression for the width Ts 



[^r]td _ £5 
2 ~ 2 



l±r] 



dS 
dE 



(63) 
(64) 

(65) 



E=[Mji]TD 



where + (— ) of ± is for the maximum(minimum) of the TD suggesting the pole on up- (pu-) 
sheet. We have found that the TD method, defined by Eq. (ITTl) and Eq. fl65l) can only find 
the shadow poles on the up-sheet which are given by exact expression Eq. fl55l) . The results 
from using the same parameters specified above are also shown in FigJH The dashed curves 
are from the exact expression Eq. (135]) and the solid dots are from applying Eq. fll7p and 
Eq. (1651) to search numerically for the poles from Eq. (135]) . Clearly TD method works very 
well in finding the shadow poles on up-sheet. In the same figure the triangles are from using 
Kelkar's prescription Eq. (1191) . Obviously, our formula Eq. (I65p works better. We have also 
examined TD for 7_ < 0. In this case the TD becomes negative and has minimum. We 
apply Eq. f|T7]) for the minimum of TD and find that TD works also well in finding the pole 
on pu-sheet. 



IV. ANALYTIC CONTINUATION OF RESONANCE MODELS 

With the analysis presented in the previous section, it is clear that the empirical partial- 
wave amplitudes determined from experimental data can not be blindly used to extract 
resonance parameters by using SP or TD methods. To make progress, one needs to construct 
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a reaction model to fit the data and then extract the resonance parameters by analytic 
continuation within the model. In this paper, we focus on a dynamical model [ol] which 
accounts for the main features of meson production reactions in the nucleon resonance region. 
Our task in this section is to develop numerical methods for finding the resonance poles from 
such models which do not have analytical forms of their solutions. We will first consider 
the simplest one-channel and one-resonance case, then two-channels and one-resonance, and 
finally two-channels and two-resonances cases. All of these models are exactly solvable such 
that their poles are known analytically and the developed numerical methods can be tested. 



A. One-channel, one-resonance 

To be specific, we consider the two-particle reactions defined by the following well known 
isobar Hamiltonian in the center of mass frame 

H = Ho + H', (66) 

with 

Ho = [E,{p) + E,{-p)] + |iVo) Mo (iVol , (67) 

H' = \9) {9\ , (68) 

where Mq is the mass parameter of a bare particle A'o which can decay into two particle 
states through the vertex interaction g in H', and Ei(p) = [m? is the energy of the 

i-th particle. The scattering operator is defined by 

t(E) = H' + H' H' , (69) 

which leads to the following Lippmann-Schwinger equation for the scattering amplitudes in 
each partial-wave 

tip ,p;E)=vip,p;E) + j^^ dq q ^ _ ^^(^) _ > (70) 

where the integration path Co will be specified later. The interaction in Eq. ( 1701) is 

v{p,p;E) = ^-^. (71) 
Eqs. (!70l) - (!7T!) leads to the following well known solution 



with 



S(E) = / dpp^- — -. (73) 
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FIG. 5: The shift of the singularity (open circle) of the propagator of the two-particle scattering 
equation Eq. (j70p as energy E moves from physical sheet to unphysical sheet. C( in (a) or Ci in 
(b) is the integration path for calculating the scattering amplitude with E on the unphysical plane. 



From the analysis in the previous section, the resonance poles can be found from t{E) = 
t{po,po;E) on the unphysical Riemann sheet defined by Impo < with pq denoting the 
on-shell momentum 



E = Jml + pl + J ml + pI (74) 



Obviously po is also the pole position of the propagator in Eq. fl70|) or Eq. (1731) . 

The physical scattering amplitude at a positive energy E can be obtained from Eq. flTOj) 
or Eq. fl72|) by setting E ^ E + ie with a positive e — > and choosing the integration 
contour Cq to be along the real-axis of p with < p < oo. From Eq. fITO]) it is clear that 
t{E) has a discontinuity on the positive real E 

Dis{t{E)) = t{E + ie) - t{E - ie) 

= 2TTip{pQ)v{po,po)t{E) , (75) 

where p{po) = Po-E'i(Po)-^2(Po) / E. Thus the t-matrix has a cut running along the real positive 
E. To find resonance poles, we need to find the solution of Eq. (I7n|) on the un-physical 
sheet with Imp < on which the pole po of the propagator moves into the lower p-plane, 
as shown in (a) of FigjSl From Eq. (|75i) . it is clear that the solution of Eq. ( |70l) . with 
the contour Co chosen to be on the real- axis < p < oo, will encounter the discontinuity 
and is not the solution on the unphysical sheet where we want to search for the resonance 



poles. It is well-known [30l. I31l . |32| . |33| that this difficulty can be overcame by deforming the 
integration path to the contour C[ shown in (a) of FiglS] By this the pole will not cross 
the cut and the integral is analytically continued from real positive E to the lower half of 
the unphysical E-sheet with Impo < 0. Obviously, the same solution can be obtained by 
choosing any contour which is below the pole position po, such as the contour Ci of (b) of 
FigEi 

With the solution of the form of Eq. (!72l) . the numerical procedure of finding resonance 
poles is to solve 

E-Mo- S(E) = 0. (76) 
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with 



S(E) = / dpp"- — y}^^ ^ , , (77) 
dpp \ (78) 

for E on the unphysical Riemann sheet defined by Impo < 0. To test this numerical 
procedure, let us consider the case that T,{E) defined by Eq. fl73l) can be calculated ana- 
lytically. Such an analytic form can be obtained by taking the non-relativistic kinematics 
Ei{p) + E2{p) = rui + 7712 +P^/(2/i) with fi = 77111712/ {mi + 7712) and a monopole form for 
the vertex function 

where is a cut-off parameter. The integration in of Eq. fl73l) can then be done 

exactly to give the following simple form 

E(i^) = ^^f^. (80) 

where po is defined by = ttii + 7712 + p1/{2n). If the imaginary part of po is positive 
(negative), it means that we choose the poles on physical (unphysical) sheet. Only the pole 
with Impo < on the unphysical sheet is called resonance, as discussed in the previous 
section. 

The resonance poles on the unphysical sheet (Impo < 0) obtained from using Eq. (180|) 
to solve Eq. (176|) are the solid curve displayed in Figj6] using the parameters: ttii = 777 = 
938.5 MeV, 7772 = w-tt = 139.6 MeV, Mq = 77ijv + + 600 MeV and cutoff momentum 
(5 = 800 MeV, for a range of coupling constant < A < 0.04. We next solve Eqs. (175]) 
and (1751) by choosing contour Ci illustrated in FigO The solutions are stable as far as the 
path is not too close to the pole. The solutions completely agree with the solid curve of the 
exact solution and hence are not displayed in Fig. [61 Here we also note that Eq. flHOl) also 
allows to calculate the poles by using the SP and TD methods. The results for the value of 
A = 0.01, 0.02, 0.03, 0.04 are the solid squares in Figj6l As expected we confirm the previous 
findings that both SP and TD work well for the single channel case. 



B. Two-channels, one-resonance 

The formula for two-channels, one-resonance case can be easily obtained by extending 
the equations in the previous subsection to include channel label i = 1,2. We thus have 

w ' T?\ I ' T?\ , f ^ ^ Vik{p',q)tkj{q,p;E) 

where Eknip) = [m\^ with rrikn denoting the mass of n-th particle in channel fc, and 

Vi, ip\ p; E) = ip') ^^^9 J (p) • (82) 
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FIG. 6: A dependence of the pole position of the one-channel, one-resonance model. The line 
represents the pole position of the t-matrix on the unphysical sheet. Solid squares are the pole 
positions extracted by using the SP and TD methods. 



Eqs. (I8I])-([82D leads to 



with 



^,{E) = I dpp' (84) 
J Co E - Ekiip) - Ek2{p) 

For the physical scattering amplitude at a positive energy E, Eq. (IHT]) is solved by setting 
E E + ie with a positive e — and choosing Cq along the real axis < p < oo. 
With Eq. (!83l) . the poles of the scattering amplitudes are defined by 

E-Mo-Si(E)-S2(^) = 0. (85) 

The poles from solving the above equation can be on one of the four Riemann sheets, pp, up, 
uu, and pu, as explained in section III. The numerical procedure for finding the resonance 
poles on MM-sheet is to solve Eq. (!85|) for E = -E'ii(poi) +-£'12(^01) = -^21(^02) +-£'22(^02) with 
Impoi < and Impo2 < 0. The integration path Cq is changed to Ci shown in (b) of Figl5] 
to calculate both self-energies Si(-E) and T,2{E) of Eq. (IMl) . Of course the contour Ci for 
the integration over the momentum for i-th channel must be below the pole poi defined by 
E — EiiipQi) — Ei2{poi) = 0. Here we note that for finding the poles on pw-sheet (wp-sheet), 
the contour Cq is replaced by Ci only for T,2{E) (Si(_E)). 
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To test numerical procedures and to further test the SP and TD methods, let us consider 
again the non-relativistic kinematics En^p) + Ei2{p) = Qi +p^/(2/ij) with 9j = mn + 
and /ij = miimi2/ {mn +mi2). This will allow us to find the exact solutions by choosing the 
monopole form factor 

We then have 

With Eq. (1871) . the poles defined by Eq. fl85l) can be found by solving algebraic equations. 
For numerical calculations, we consider a case similar to vrA^ scattering in 5*11 partial wave: 
(l)channel-l is vrA^ with mn = = 139.6 MeV , mi2 = rri]^ = 938.5 MeV, and f3i = 
800 MeV, (2)channel-2 is rjN with mai = = 547.45 MeV , m22 = niN = 938.5 MeV, 
and (32 = 800 MeV, (3) bare mass Mo = m.,^ + rriN- + 600 MeV. The results are shown in 
FiglTl The dash-dotted (solid) curves are the calculated poles on ww-sheet (wp-sheet) with 
the coupling constants Ai = 0.02 for No tcN and a range of A2 = — 0.02 for A^^o rjN . 
We see that when A2 is 0, which is the single-channel case, the ww-pole and the wp-pole are 
on the same position. They then split as A2 increases. 

We next evaluate Eq. (18^ for E on the ww-sheet, wp-sheet or p-u-sheet by appropriately 
choosing the path Cq, as described above. The poles are found when the calculated 
and T,2{E) satisfy Eq. (ISBj) . We find that the poles obtained by this numerical procedure 
reproduce accurately the dash-dotted (wM-sheet), solid (wp-sheet) and dashed (p-u-sheet) curve 
in FigJT] and hence are omitted there. Thus this analytic continuation method can be used 
in practice to find the resonance poles (i.e. poles on uu-sheet as defined in section III) for 
general case that T,i{E) can not be integrated analytically. 

We now turn to examining the SP and TD methods. In Figd we show that the SP 
method (open squares) reproduces the poles (dash-dotted curve) on the ww-sheet only at 
A2 < 0.015. When A2 is larger than 0.015 where the magnitude of hxvEji continues to 
increase, the speed has no maximum and the SP method fails to find the pole. This is 
another example showing that SP method has its limitation. On the other hand, the TD 
method (solid squares) can reproduce the poles on the wp-sheet (solid curve) and pw-sheet 
(dashed curve) in the considered range of parameters. Here we see that the SP and TD find 
different pol es which have different physical meanings in the Hamiltonian formulation. One 
can show[l5l| that the poles on ww-sheet are due to the process that an unstable system is 
created and then decays during the collision and are called the resonance poles. The physical 
interpretations of the poles on pu and -up-sheets remain to be developed. 

It is interesting to point out here that the two-channel Briet-Wigner form analyzed in 
detail in the previous section can be derived from the two-channels, one-resonance model if 
the non-relativistic kinematics is used. To see this, we first write the non-relativistic relation 
between the S-matrix and the T-matrix 

Sij{E) = 6ij + 2iTij{E), (88) 
Tij{E) = -rf^yniPiHjPjtij{pi,pj; E), (89) 

where Pi is the on-shell momentum in channel i 



p, = j2fii{E-Qi). (90) 
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FIG. 7: A2 dependence of the pole positions on the itn-sheet (dash-dotted curve), up-sheet (sohd 
curve) and pn-sheet (dashed curve) of the two-channels, one-resonance model. The open squares 
(solid squares) are obtained from using the SP (TD) method. 



With the above and the analytic form Eq. (1871) for Sj(£') , we can write the 1 — >■ 1 elastic 
scattering amplitude of Eq. fl89l) as 

Tu{pi,Pi,E) = . -— , (91) 

E - M{E) + zpi7i(pi) + tp2l2{P2> 



where ji^pi) = 7ifiigf{pi) > and 

M{E) = M' + Y.P I P'dP^ ^ (92) 

k J E - ruki - mk2 - ^ 

where P means taking the principal- value integration. By using Eq. fl88l) . we then have the 
1^1 elastic part of the S-matrix 

g ^ E - M{E) - ipaijpi) + ip2l2{p2) .ggx 
11 E-M{E)+tp^^M+mi2{P2)' ^ ^ 

If the E-dependence of M{E) and jii^Pi) are further neglected, Eqs. (I9T1) and (!93l) are identi- 
cal to what are usually called the two-channel Breit-Wigner resonant amplitude discussed in 
the previous section. Thus the conditions under which SP and TD are valid can be related 
now to the parameters of the vertex function gi{p) within this two-channels, one-resonance 
model. 
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C. Two-channels, two-resonances 



For the two-channels, two-resonances case, the scattering amplitude is defined by the 
same Eq. flHTl) . but with the following driving term 

1 1 

Vij{p\p] E) = 9nip') ^ _ j^^ 9ji{p) + 9i2{p) ^ _ ^^ 9j2{p). (94) 

The scattering amplitude is then of the following form 

ti,{p',p; E) =Y,g,,^{p')[D-\E)]^,pgjM- (95) 

The propagator D{E) in Eq. fl95|l is 

[D{E)]^^p =[E- M^%,f3 - ^aAE), (96) 

with 

'^^ ' ^Jco E-Eaiq)-EAq)+te ^ ^ 

The poles are defined by 

BetDiE) = [E-Mi-Sn(E)][E-M2-S22(^)] -Si2(^)S2i(^) 

= 0. (98) 

The numerical procedures of finding the resonance poles on the wu-sheet from Eqs. f PT]) 
and fl^S]) are the same as that in the previous two subsections. Namely the path Co of Eq. 
flU7|) is set to be the path Ci shown in Fig IS] in evaluating the integrals for E on the uu-sheet 
where the on-shell momenta are Impj < for i = 1, 2 channels. To test this, we again choose 
the non-relativistic kinematics and the monopole form factor like Eq. (1861) . The self energy 
Tia,i3{E) then takes the analytic form similar to Eq. (IHOl) . The condition Eq. (p8|) can then 
be expressed in a analytic form from which the pole positions on the unphysical sheet can 
be easily obtained. 

We only state that the resulting poles on the unphysical sheets are reproduced by the 
numerical analytic continuation method described above. Instead our focus here is to further 
test the SP and TD methods for the situation that two resonances are close and could overlap. 
We again consider irN and riN channels and use the following form factor 

9^M = YVM^ ^''^ 



where a = 1,2 denote the a-th bare state with mass Mq. The four cut-off parameter 
and four coupling constants \ia are taken to be: Pn = = P21 = P22 = 800 MeV, 
All = 0.005, A12 = 0.01, A21 = 0.003 and A22 = 0.008. One of the bare masses Mi is fixed in 
the calculations. M2 is defined by 

M2 = M^ + Mn + M2, (100) 

where M2 is varied for examining how the poles move as M2 moves away from the Mi = 
M^ + Mn + 550 = 1628 MeV. 
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The pole positions are searched numerically by using the analytic continuation method 
described above. There are two poles on the uu-sheet and the other two on the up-sheet. 
As M2 varies, these two poles will develop two trajectories. They are the crosses connected 
by the solid curves shown in FigJS]for uu-sheet and FiglHlfor the up-sheet. According to the 
findings we made in section III, the poles, {Mji, —iTji/2), found by the SP (TD) method 
should be compared with the poles on uu (up) sheet of Fig. [S] ©• We now discuss these 
two comparisons. 

We see from Fig. [8] that in the regions near M2 = 700 MeV, the positions (ReE) of these 
two poles are far from each other and we find that SP_ (open squares connected by dashed 
lines near the point marked 700) works well. When M2 is reduced to 600 MeV where the 
positions (ReE) of two poles move closer, the SP method can find only one pole ( open 
squares connected by dashed line) near the top end of the trajectory on the right hand side. 
Apart from the points on the dashed lines, SP method fails to find poles close to the poles 
on the solid curves which are obtained numerically by the analytic continuation method. 

The results for examining TD is shown in FiglHl We see that TD can find two poles on 
the up-sheet in the considered range of M2 = 600 — 700. The results are the open squares 
connected by dashed lines which are indistinguishable from the crosses connected by solid 
curves which were obtained by analytic continuation method. But TD obtains another two 
poles at M2 = 600 and M2 = 625, as indicated by the dashed line in the middle of Figj9l 
The positions of these two poles are very close to those obtained from SP and they are 
interpreted as poles on uu-sheet. As we have discussed in section II, TD is sensitive to both 
zero and pole of the S-matrix on pp and uu-sheets. In this example, width of the poles on 
MM-sheet become comparable to the poles on Mp-sheet(zero on physical sheet) and hence TD 
could find pole on uu-sheet. The results shown in Fig. |S] and Fig. [H] further indicate the 
limitation of SP and TD methods. 

The results from the above several models have shown that the TD method based on Eq. 
f|T6|) . where the phase of the elastic channel is used instead of the eigen phases discussed 



in Refs. [23|, |24, l25j, gives both the resonance poles on the uu-sheet and the zeros(shadow 
poles). Our findings could provide some information for investigating the differences between 
Refs. [25] andfsj]. We also want to mention here that an improved SP method using higher 
order derivatives of the amplitudes was proposed in Ref. |35i]. It may be interesting to 
compare this method with the TD method. While it could be interesting to address the 
questions concerning these recent developments, they are far from the main focus of this 
paper and will not be discussed further. 



V. ANALYTIC CONTINUATION OF RESONANCE MODEL WITH UNSTABLE 
PARTICLE CHANNELS 

For meson-baryon reactions, the nucleon resonances can decay into some unstable particle 
channels such as the vrA, pN, aN considered in the model of Ref.j^. Here we discuss the 
analytic continuation method to find resonance poles within such a reaction model. 

It is sufficient to consider the one-channel and one resonance case. The scattering formula 
is then identical to that presented in subsection IV.A. The only difference is that one of the 
particles in the open channel can further decay into a two particle state. To be specific, 
let us consider the vrA channel. Within the same Hamiltonian formulation [91 used in the 
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FIG. 8: Pole positions (crosses connected by solid lines) on the un-sheet of the two-channels, 
two-resonances model. The results from using the SP method are the open squares connected by 
the dashed lines. The numbers on the figure are the value of M^. 
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FIG. 9: Pole positions (crosses connected by solid lines) on the up-sheet of the two-channels, 
two-resonances model. The results from using the TD method are the open squares connected by 
the dashed lines. The numbers on the figure are the value of M2. 
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previous section, the scattering amplitude can then be written as 



t{p',p,E) 



with 



^nA{E) = p'^dp 

JC2 



gN*,7TA{p')gN*,7TA{p) 

E-Mo- ^^E) 



where 



^a{p,E) 



q^dq 



E-E^{p)-Ea{p)-^a{p,E) 
E - E^{p) - [{EM + EN{q)Y+p^Y'^ ■ 



(101) 



(102) 



(103) 



To obtain the vrA self energy for complex E, the analytic structure of the integrand of Eq. 
(11021) should be examined first. The discontinuity of the vrA propagator in the integrand of 
Eq. (11021) is the vrvrA^ cut along the real axis between ±po (— Po < P < po) which is obtained 
by solving 



E = E^ipo) + [(m, + niMf + Pq]'^"- 



(104) 



For finding the resonance poles on wu-sheet with Impo < 0, the integration contour C2 of 
Eq. (11021) must be chosen below this cut which is the dashed line in FigJTOl There is also a 
singularity in the integrand of Eq. (I102p at momentum p = Px, which satisfies 



E - EM - Ea{p.) - Sa(Px, E) = 0. 



(105) 



Physically, this singularity corresponds to the vrA two-body scattering state. For E with 
large imaginary part, px can be below the ttttN cut as also indicated in FigJTOl Therefore 
the integration contour of momentum p must be chosen to be below the TryriV cut (dashed 
line) and the singularity p^, such as the contour C2 shown in FigJTUl 

The singularity position go of the propagator in Eq. (I103p depends on spectator momen- 
tum p 



E - E^ = [(E.(go) + ^iv(go))' + P 



211/2 



(106) 



Therefore the singularity go moves along the dashed curve in FigHT] when the momentum p 
varies along the path C2 of FigJTUl To analytically continue Tia{p,E) from positive energy 
E to the un-physical plane with Imp < 0, we need to choose the contour C3 of Eq. (I103P 
which must be below go. A possible contour C3 is the solid curve in FigJTTl 

To verify the numerical procedures described above, we again consider non-relativistic 
kinematics and monopole form factor. With the similar analytic form Eq. (IHOl) . we have 



^a{p,E) 



2 aS 

'^f^7TN9A,TTNyA,TTN 

'^{k + lPA,.N? 



(107) 



where 



2/i^7v E - 2m^ 



rriN 



p^ 



.1/2 



2/i 



TTTvN , 



(108) 
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FIG. 10: Contour C2 for calculating the vrA self energy on unphysical sheet. See the text for the 
explanations of the dashed line and the singularity px- 
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FIG. 11: Contour C3 for calculating the vrA^ self energy on the unphysical sheet. Dashed curve is 
the singularity qo of the propagator in Eq. (jl03p . which depends on the spectator momentum p on 
the contour C2 of FiglTUl 



with ^ttttN = fTiT^imT, + mAr)/(2m^ + mjq). With Eq. (11071) . we can solve Eq. (11051) and 
verify its relation with vrvrA^ cut as discussed above and illustrated in FigHDl Eq. (11071) and 
the chosen monopole form factor also allow us to get 



with 



2(/c + iPA,nNr 
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FIG. 12: N* Green function. The solid (dash-dotted) curve is the real (imaginary) part of exact 
Green function Gjv*(-E) = 1/{E — Mq — Y!,t^/\{E)). They are compared with the dashed (dotted) 
curve of the real (imaginary) part of the approximate Green function G°^°^ ■ 



Unfortunately, Eq. ( 1109P can not be integrated out analytically for directly checking our 
numerical procedure for searching resonance poles. 

We test our analytic continuation method by the following procedure. We calculate Eq. 
(11091) numerically to find the pole position Eji by solving Eji — Mq — S^aI-^/j) = of the 
denominator of Eq. (llOip . With the parameters: 

pN',nA = SOOMeV , gN',nA = 0.02MeF~^/^ 
pA,nN = 200MeV , gA,nN = 0.05Me\/-^/^ 
Mo = Q50MeV + m^ + mN , 

we find Er = (1679.1, — 33.6z) MeV. We then construct an approximate propagator 

Gr™^(i?) = . (Ill) 

For the positive E, we find that G'^f^"'^ (E) is in good agreement with the direct calculation 
of Gn*{E) = l/{E-Mo- ^nA{E)) by using Eq. ffT09D . The resuhs are shown in FigUl It 
is clear that the resonance pole found by our analytic continuation method can reproduce 
what is expected for a resonance propagator for real positive E. In this way our numerical 
procedure is justified and can be applied to solve Eqs. (110 II) - (11031) . 



VI. RESONANCE MODEL WITH NON-RESONANT INTERACTIONS 

With the numerical methods described above, we can proceed to extract the resonance 
poles within a coupled-channels model which also include non-resonant interactions. In this 
section, we explain how this can be done for the irN model developed in Refs.jsl, [2^ . 

Recalling the formulations presented in Refs.jsl, [22|, the t-matrix considered is of the 
following form 

tijip',p; W) = u,{p\p- W) + Ef,.(p'; W)\D-\w%,pf,^p{p- W). (112) 
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where i,j can be stable channels vrA^ and rjN, or unstable channels irA^pN, and aN, and 
a denotes a bare resonant state with a mass Ma- 

For extracting resonance poles, we now apply the methods presented in previous sections 
to choose appropriate contours for calculating various integrations on unphysical E-plane 
with Imp < 0. The non-resonant t-matrix tij{p',p; E) is defined by the following couple- 
channel equations with the non- resonant potential Vij{p',p), 

r / / { ' \ , f ^ 2 Vikip' , q)ikj{q, P] E) 



k 



where the contour is C4 = Ci ((b) of Figj5]) for k = irN^rfN, and C4 = C2 (FigiTO]) for 
k = 7^A,pN,aN, and 



Ek{p) = \/ml, + p2 + ^ml^ +p2_ (114) 

The self-energies in Eq. flllSp are S^Ar(g, E) = T,j^^{q, E) = and T,T,^{q, E), T,p^{p^ E), and 
T,„N{q,E) are defined by the same Eq. ( 1102^ with appropriate changes of mass parameters 
and the choice of contour C2 and C3 shown in Figs. [TOlandfTTl 

The dressed vertex Ti{p] W) is determined by the bare vertex Ti{p) and the meson-baryon 
loop, 

r f r { \ I ^ f ^ 2 iik{p,q;E)Tk{q;E) 

Ti.a p; E) = Ti^aiP) + 2^ dqq — , • • 

^Jci E - Ek{q) -T.k{p,E) + ie 

The resonance propagator D{E) in Eq. flll2p is 

[D{E)]a,p =[E- Ma]5a,p - ^aAE), (116) 

with 

^aAE) = 2^ dqq — fa , • • 

^JCi E - Ek{q) - T.k{q, E) + te 

In Ref. 23], the above equations are solved on the real-axis by using the standard method 
of subtraction. Here we solve the equations by choosing contours indicated above. We first 



verify that our numerical results obtained here for positive real E agree with that of Ref. [22 
This establish our numerical procedure in this complex five-channel model. 

Here we show the results for some of the S, P, D, and E partial waves of vrA^ scatter- 
ing within the model of Ref.[22|. We search the resonance poles by looking for zeros of 
the resonant propagator D{E) defined by Eq. flll6p . Our results from using the analytic 
continuation method are shown in the second column of Table [U They are compared with 
those extracted by using SP and TD methods described in the previous sections as well as 
the values listed by PDG. The Breit-Wigner resonance parameters are also given by PDG, 
but are not considered here. We see in Table [T] that the SP method fails to find the first 
resonance pole in the Su partial wave. We also see that while the real parts of the resonance 
poles from different approaches are within the ranges of PDG, the extracted imaginary parts 
can differ by as much as a factor of two or three. 

The model of Ref. [2^] is currently being improved by also fitting other data of tcN and tA^ 
reactions. For example, some progress has been made to also fit the data of ttN 7r7riV[36|, 
7A^ 71 N 3^, and ttN — »• ?7iv[3^. We thus do not include the results for other partial 
waves in Table [B Our purpose here is to simply demonstrate how the analytic continuation 
method works for a realistic model. The full resonance parameters, including the extracted 
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TABLE I: Resonance poles extracted from the ttN scattering amplitudes of Ref. 





Analytic Continuation 


Speed PLot 


Time Delay 


PDG 




Re 


Im 


Re 


Im 


Re 


Im 


Re 


Im 


Sll 


1540 


-191 






1543 


-52 


1490 ~ 1530 


-45 ~ -125 




1642 


-41 


1644 


-89 


1645 


-61 


1640 ~ 1670 


-75 ~ -90 


S31 


1563 


-95 


1574 


-67 


1616 


-53 


1590 ~ 1610 


-57 ~-60 


P33 


1211 


-50 


1212 


-49 


1212 


-49 


1209 ~ 1211 


-49 ~-51 


D13 


1521 


-58 


1525 


-57 


1522 


-11 


1505 ~ 1515 


-52 ~-60 


F15 


1674 


-53 


1671 


-59 


1683 


-24 


1665 ~ 1680 


-55 ~-68 



residues and the relations to the Breit-Wigner parameters Usted by PDG, extracted from 
our complete analysis of all irN'j, N ttN, rjN, ttttN will be reported elasewhere. 

We now turn to discussing whether the extracted resonance pole M can be used to 
evaluate the A^* propagator defined by Eq. flll6p for the physical positive E. Let us consider 
the 5*31 case listed in Table HI Its A^* propagator Gn*{E) can be written as 



E(E) 



1 



E-Mo-S(E)' 



(118) 
(119) 



By using the analytic continuation methods described in the previous sections, the resonance 
energy M = (1563 — i95) MeV is found numerically by solving 

M-Mq- E(M) = 0. (120) 

We now perform the Laurent expansion of Gn*{E) for real E around the pole position 



M 



GN'iE) = [(1 



S'(M))(^ - M) - is"(M)(E - Mf 



E-M 
1 

E-M 
1 



1 - S'(M) - iS"(M)(E - M) + 
1 ^S"(M)(E-M)^ 



l-E'(M) 2(1-S'(E))2 
1 S"(M) 
E-M ' 1 - S'(M) ^ 2(1 -E'(M))2 



+ 



;i21) 



The naive 1/{E — M) works well for a model studied in the previous section. However 
when the enery dependence of the self energy S becomes large, two important modifications 
should be considered: (1) at the pole, the residue is not one but modified by the field 
renormalization factor Z = 1/(1 — S'(M)). (2)The second term of the last expression gives 
a constant term. 

With the above expansion Eq. (11211) . we can introduce three different approximate forms 
for the N* propagator 



E-M' 



;i22) 
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FIG. 13: Comparisons of various resonance propagators: exact propagator Gn*{E) (cross), 
G^^l{E) (dashed curve ) , G^^l{E) (dash-dotted curve) and G^^l{E) (soUd curve). The real (imagi- 
nary) parts are shown in the (a) ((b)) parts of the figure. 



= E^-l-S'(M)' ^^^^^ 

G^'^ (E) = ^— ^ + ^^^(^) n 24) 

' E-M l-S'(M) 2(1-S'(M))2- ^ ' 

In FigIT3| we compare the above three approximate propagators with the exact result of 
Gn4E) of Eq. ffTTSD . The simple gP,{E) (dashed curves) of Eq. ffT22D are far from the 
exact green function Gn*{E) (cross) defined by Eq. flllSp . When the factor 1/(1 — E'(M)) is 
included, we obtain the dashed-dotted curves forGS^l(E). The solid curves are from G^^l{E). 
We see that the constant term 2(i-s'''|m))2 -'^l- P24l) mainly affects the imaginary part of 
the propagator. 

Eq. (11211) shows that the simple pole approximation G^^l{E) = ^^j , which can be cast 
into the usual Breit-Wigner form 1/{E - MR + irR/2) with Mr = Re M and Tr = -2ImM, 
works poorly. We also find that the pole parametrization of the resonant propagator using 
the poles extracted by the TD and SP methods also work poorly. For the considered 5*31 
case, the pole positions from using these two methods are: Esp = (1574, —67i) and Etd = 
(1616, —53?), as given in Tabled In Fig{Tll we compare the exact propagator (cross) Gn* (E) 
of Eq. (11181) with the following two propagators 

G^^ME) = ^ , (125) 

G^^ME) = ^ • (126) 

Clearly, phenomeno logical forms Eqs. fll25p - (11261) can not account for the complex coupled- 
channel resonant mechanisms. 
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FIG. 14: The resonant propagator Gn*{E) (cross) is compared with Gjy* ^sp{E) (sohd curve) using 
the SP poles and Gn* ,td{E) (dashed curve) using the TD poles. The real (imaginary) parts are 
shown in the (a) ((b)) parts of the figure. 



VII. SUMMARY 

In this paper, we have presented a pedagogical study of the commonly used Speed- 
Plot (SP) and Time-Delay (TD) methods for extracting the resonance parameters from the 
empirically determined partial-wave amplitudes. Using a two-channel Breit-Wigner form 
of the S-matrix, we show that the poles extracted by using theses two methods are on 
different Riemann sheets. The SP method can find resonance poles on the unphysical uu- 
sheet, while the TD method can find poles and zeros of S-matrix on uu or pp-sheets and 
therefore its validity is sensitive to the poles on up or p-u-sheets. Furthermore, we also show 
numerically that these two methods can fail to find those poles. Our results support the 
previous findings that these two methods must be used with cautions in searching for nucleon 
resonances from the meson-nucleon reaction data in the region where the coupled- channel 
effects are important. 

We then develop an analytic continuation method for extracting the resonance poles 
within a Hamiltonian formulation of meson-nucleon reactions. The main focus is on resolving 
the complications due to the coupling with the unstable vrA, pN , and aN channels which 
can decay into t^ttN states. Explicit numerical procedures are presented and verified within 
several exactly solvable models. The results from these models are also used to further 
demonstrate the limitation of the SP and TD methods. 

As a first application of the developed analytic continuation method, we present the 
results from analyzing the ^n, 5*31, P33, D13 and F}^ amplitudes of the dynamical coupled- 
channels model of vrA^ reactions developed in Ref . [2^ . We also analyze the resonance prop- 
agators and show that the simple pole parametrization of the resonant propagator using the 
poles extracted from SP and TD methods works poorly. 

With the progress made in this work, we can proceed to extract all nucleon resonance 
parameters within the model of Ref. 2^. However, this can be done more accurately only 



when the coupling with the unstable vrA, pN and aN channels are better determined by 
also fitting the two-pion production data. Our progress in this direction will be reported 
elsewhere. 
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